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We use the framework of a relativistic constituent quark model to study the semileptonic 
transitions of the Be meson into (cc) charmonium states where (cc) — r)c CSo), J/tp CSi), 
XcoCPo), Xci (^A), hcCPi), XC2CP2), Tp{^D2). We compute the g^-dependence of aU rele- 
vant form factors and give predictions for their semileptonic Be decay modes including also their 
r-modes. We derive a formula for the polar angle distribution of the charged lepton in the (liyi) cm. 
frame and compute the partial helicity rates that multiply the angular factors in the decay distri- 
bution. For the discovery channel Be J/'il){-^ fi^ ijl~)Iv we compute the transverse/longitudinal 
composition of the J/?/) which can be determined by an angular analysis of the decay J/ip 
We compare our results with the results of other calculations. 
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I. INTRODUCTION 

In 1998 the CDF Collaboration reported on the observation of the bottom-charm Be meson at Fermilab [1]. The 
Be mesons were found in an analysis of the semileptonic decays Be J/ipli' with the J/^p decaying into muon 
pairs. Values for the mass and the lifetime of the Be meson were given as M{Be) — 6.40 ± 0.39 ± 0.13 GeV and 
T{Be) = 0.46^0 i6(stat)±0.03(syst) • 10""'^^ s, respectively. First Be mesons are now starting to be seen also in the Run 
II data from the Tevatron [2, 3]. Much larger samples of Be mesons and more information on their decay properties are 
expected from the current Run II at the Tevatron and future experiments at the LHC starting in 2007. In particular 
this holds true for the dedicated detectors BTeV and LHCB which are especially designed for the analysis of B physics 
where one expects to see up to 10^° Be events per year. 

The study of the Be meson is of great interest due to some of its outstanding features. It is the lowest bound state 
of two heavy quarks (charm and bottom) with open (explicit) flavor. As far as the bound state characteristics are 
concerned the Be meson is quite similar to the J^'^ = states rje and r]t in the charmonium (cc-bound state) and 
the bottomium (66-bound state) sector. However, the r]e and rjh have hidden (implicit) flavor and decay strongly and 
electromagnetically whereas the Bc-nieson decays weakly since it lies below the i3Z)-threshold. 

The Be meson and its decays have been widely studied in the literature. The theoretical status of the Sc-meson was 
reviewed in [4] . The Be lifetime and decays were studied in the pioneering paper [5] . The exclusive semileptonic and 
nonleptonic (assuming factorization) decays of the Bc-nicson were calculated in a potential model approach [6] . The 
binding energy and the wave function of the Sc-meson were computed by using a flavor-independent potential with 
the parameters fixed by the cc and bb spectra and decays. The same processes were also studied in the framework 
of the Bethe-Salpeter equation in [7] , and, in the relativistic constituent quark model formulated on the light-front 
in [8]. Three-point sum rules of QCD and NRQCD were analyzed in [9] and [10] to obtain the form factors of the 
semileptonic decays of B'^ J/ip{r]e)l~^v and B^ Bs{Bl)l'^v. As shown by the authors of [11], the form factors 
parameterizing the Be semileptonic matrix elements can be related to a smaller set of form factors if one exploits the 
decoupling of the spin of the heavy quarks in the Be and in the mesons produced in the semileptonic decays. The 
reduced form factors can be evaluated as an overlap integral of the meson wave- functions which can be obtained, for 
example, using a relativistic potential model. This was done in [12], where the Be semileptonic form factors were 
computed and predictions for semileptonic and non-leptonic decay modes were given. 

In [13] we focused on its exclusive leptonic and semileptonic decays which are sensitive to the description of long 
distance effects. From the semileptonic decays one can obtain results on the corresponding two-body non-leptonic 
decay processes in the so-called factorization approximation. The calculations have been done within our relativistic 
constituent quark model based on an effective Lagrangian describing the coupling of hadrons H to their constituent 
quarks. The relevant coupling strength is determined by the compositeness condition Zh = [14, 15] where Zh is 
the wave function renormalization constant of the hadron H . 

The relativistic constituent quark model was also employed in a calculation of the exclusive rare decays Be — > 
D{D*)ll [16] and of the nonleptonic decays Be DJW and Be 0^0° [18]. In the latter case we confirmed that 
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the nonleptonic decays Be — > DgD^ and Be DgD^ are well suited to extract the CKM angle 7 through amphtude 
relations, as was originally proposed in [19, 20]. The reason is that the branching fractions into the two channels are 
of the same order of magnitude. 

In this paper we continue the study of Be decay properties and calculate the branching rates of the semileptonic 
decays ^ (cc)?!/ with (cc) =r;e(i5o), J/V'C^i), Xco (^^o), Xci C'A), Ki^Pi), Xc2 ('^2), (^Da). We compare 
our results with the results of [6, 21] where it was shown that these decay rates are quite sizable and may be accessible 
in RUN II of Tevatron and/or the LHC. Two-particle decays of the i?c-meson into charmonium states have been 
studied before in [22] by using the factorization of hard and soft contributions. The weak decays of the Bc-meson to 
charmonium have been studied in the framework of the relativistic quark model based on the quasipotcntial approach 
in [23]. In this paper we compute all form factors of the above semileptonic i?c-transitions and give predictions 
for various semileptonic B^ decay modes including their r-modes. Prom a general point of view we would like to 
remark that the semileptonic decays of the r-lcpton have been studied within perturbative QCD. It has allowed one 
to determine the strong coupling constant with a high accuracy (see e.g. [24]). We have improved on our previous 
calculation [13] in that we no longer employ the so-called impulse approximation. In the impulse approximation one 
assumes that the vertex functions depend only on the loop momentum flowing through the vertex. Dropping the 
impulse approximation means that the vertex function can also depend on outer momenta according to the flow of 
momentum through the vertex. A comparison with the results for the decays into the para- and ortho-charmonium 
states (cc) = r]c (^Sq), J/il) i^Si) [13], which was done in the impulse approximation, shows a « 10% downward effect 
in the rates when the impulse approximation is dropped. 

II. BOUND STATE REPRESENTATION OF THE CHARMONIUM STATES 

The charmonium states treated in this paper are listed in Table I. We have also included the purported I^-wave 
state V(3836) whose quantum numbers have not been established yet. Table I also contains the quark currents used 
to describe the coupling of the respective charmonium states to the charm quarks. The masses of the charmonium 
states listed in Table I are taken from the PDG [25]. 



TABLE I: The charmonium states '^^'^^L j. We use the notation d=d — d- 



quantum number 


name 


quark current 


mass (GeV) 


jPC ^ Q-+ (s^Q^ L^Q) 


'So = Vc 




2.980 


J^'^- = i (S=i, L=0) 


■'.S'i = J/r 


q-'' q 


:!.()!)r 


jPC ^ Q + + (g^;^^ L^;^) 


^Po = XcO 


qq 


3.415 


jPC ^ (S^;^^ L = l) 




qy^y^q 


3.511 


J^^ = 1+- (S=0, L=l) 


Ipl = h,{lp) 


q d 7 <? 


3.526 


jPC ^2++ (S=l, L=l) 


= Xc2 


{i/2)q (^7'' 9" -^7" a") q 


3.557 


J^c ^ 2~ (S=l, L=2) 


^Di = V(3836) 


(i/2)g(7V d'^ +7V a") g 


3.836 



Next we write down the Lagrangian describing the interaction of the charmonium fields with the quark currents. 
We give also the definition of the one- loop self-energy or mass insertions (called mass functions in the following) 11 (p^) 
of the relevant charmonium fields. 

We can be quite brief in the presentation of the technical details of our calculation since it is patterned after the 
calculation presented in [13] which contains more calculational details. We treat the different spin cases {S = 0, 1,2) 
in turn. 



Spin S=0: 



U{x - y) 



= l^{x){^-m'')^{x)+g(f){x)Jg{x), 
= tg^T{Jg{x)Jg{y)})o, 



(1) 



Z = 1- n'(m2) = 1 - ^n^(m2) = 0, 
Jqi^) = jjdxidx2Fcc{x,Xi,X2)q{xi)Tq{x2) (r = i7^,/), 

Fcc{x, Xi,X2) = 6 ^X - ^^^ ^ $cc {{Xl - X2f) , 

n'(m^) is the derivative of the mass function n(p^). 



Spin S=l: 



1 

2 

dij,(j/^{x) = (leaving three independent components), 



Cs=i{x) = -^-<i>^,{x){U^w?W{x)+g<i>^,{x)J'^^{x), 



W^^{x-y) = -i5'(T{JA'(a;)J,-(2/)})o, 

W'^ip) = J d'^xe-'P''W'{x)=g^"'n'^^\p'^)+p^'p''U^'^\p'^), 



•^9(2;) = j j dxidx2Fcc{x,xi,X2)q{xi)T^ q{x2), 

The spin 1 polarization vector e^^^ {p) satisfies the constraints: 

e^^^ (p) p** = transversality, 

Yl ^f^Pyl^^^P) = -dill' + ^^^ completeness, 
A=o,± ™ 

ej^^'^'^e^'^ = -^AA' orthonormality. 



Spin S=2: 

£s=2(a;) = \ cl,^,{x){U - m^)r"{x) +9^^.{x) Ji^"{x). 

(a;) = 4'"'^{x), dij,4i'^"{x) = Q, (l)^{x) = Q, (leaving 5 independent components), 

n''-.°/3(x-y) = <T{Ji^''{x)jf{y)]>o, 

n^.,a/3(p) = Jd^xe-'P''W'''^{x) = ^{g^"g-''^+g^'^g'"') U^^^ip^) 
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J^"{x) = J|dxldx2F,,{x,XuX2)q{xl)T^'''q{x2), 

r^" = \ T +7" d'^) , ^ (7V +7V a") . 

The spin 2 polarization vector e^'^'(p) satisfies the constraints: 



j2 ^'■^^^''^-^^ 



A=0,±l,±2 
t(A) (A')m<^ = 



= Sw 



symmetry, 
transversality, 
tracelessness, 
completeness, 

orthonormality. 



where 



We use the the local representation for the quark propagator when calculating the Fourier-tansforms of the mass 
functions. The local quark propagator is given by 



S,{x-y) = {T{q{x)q{y)})o = J e'^^^"^) S,{k) = ^ 

For the mass functions one needs to calculate the integrals 



(4) 



no(p') 

T{P,S) 



~ I ^i^lci-k')Tr[TS{k-p/2)rS{k+p/2) 
^7^ /. 



r dH 
~ J ^ 



^i^{-k^)Tr V S{k-p/2)T'' S{k+p/2) 



nr(p) 

T^'{V,A,PV) = 7'^, 7V, 2ifcV- 

fi-^'^^ip) = I ^cc(-fc') Tr fP"'' S{k - p/2) S{k+p/2) 
J in I L 

Yt^^{T, PT) = i (7" k"" + i' k^) , i {j^-f^ k" + Y-i^ k") . 

The functional form of the vertex function $cc(^fc^) ^ind the quark propagators Sg{k) can in principle be determined 
from an analysis of the Bethe-Salpeter and Dyson- Schwingcr equations as was done e.g. in [26]. In this paper, 
however, we choose a phenomenological approach where the vertex functions are modelled by a Gaussian form, the 
size parameters of which are determined by a fit to the leptonic and radiative decays of the lowest lying charm and 
bottom mesons. For the quark propagators we use the above local representation Eq. (4). 

We represent the vertex function by 



^cc{-k^) 



„Scc fe 



1 
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where Acc parametrizes the size of the charmonium state. The quark propagator can easily be calculated using 
Feynman parametrization. One has 

oo 

Sgik±p/2) = {mq+f6±ti/2) j dae-°'^<-^^-P''^^^K 



We then transform to new a-variables according to 
and make use of the identity 

oo oo oo 

j j daida2 f {011,012) = J dtt JJ daida2S {1 — ai — 02) f{tai,ta2) ■ 



One then obtains 



where 







Ct 



t 



2(l+t)Scc' " 1 + t 

The symbol < ... > stands for the two-fold integral 

00 1 



a-- 



(•••) = dt 



t 







da e 



(•••), 



z = t [m^ — a(l — a)p^] 



(5) 



It is then convenient to shift the loop momentum according to —> — Cp p. The ensuing tensor integrals can be 
expressed by scalar integrals according to 

[ ^fi-k')k'^k-^ = -,9'''' I ^f{-k')k\ 



The remaining scalar integrals can be integrated to give 

J4 



J TT t 



(6) 



For the mass functions one finally obtains 



97^7^ q 
q d Tq 
q^ (d'^Y+Tl" 



U{p^)p = {2ct + ml + {l/A-cl)p^) 

n{p^)s = {2ct - ml + a/A - cl)p^} 
U{p^)v = {ct+ml + {l/A-cl)p') 
n{p^)A = {ct~ml + {l/A-cl)p^) 
n{p'')pv = 2ct{Sct+ml + (1/4 - c2) p2 ) 
II{p^)t = 2ct{ct + ml + (1/4 - c^) p^ ) 
Uip')pT = 2ct{ct-ml + (1/4 - cl)p^ ) 
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The mass functions 11/ (p^) enter the compositeness condition in the derivative form 

where the prime denotes differentiation with respect to p^. The differentiation of the mass functions result in 

n'(p')p = (-2sec-S[2ct + m2 + (l/4-c2)p2] + i/4-c2) 

n'(p')s = (-2sec5[2ct-m2 + (l/4-c^)p2] + i/4_c2 ) 

n'(p'V = (-2scc5[ct + m^ + (l/4-c2)p2] + i/4_c2 ) 

n'(p')A = (-2sec5[ct-m^ + (l/4-c^)p2] + i/4-c2) 

n'(p')py = 2ct{-2s^z[^Ct + ml + {l/A-cl)p^] + l/A-cl) 

n'(p2)y = 2ct(-2sec^[ct + m2 + (l/4-c^)p2] + i/4_c2 ) 

n'(p2)pT = 2ct(-2sec5[ct-m2 + (l/4-c^)p2] + i/4_c2 ) 



where 

z = —ta(l — a) — -— — I Q — 
^ ' l+t \ 2 

III. THE SEMILEPTONIC DECAYS Be (CC) +L + P 

Let us first write down the interaction Lagrangian which we need for the calculation of the matrix elements of the 
semileptonic decays Be {cc) + 1 + P. One has 



Ant (a;) = 56c {x) ■ J^+(x) + C/cc (Pcc{x) ■ Jcc{x) + Vbc {cO'' b) ■ {I v), 



Gp 

Jbci^) = JJ dxidx2F\,c{x,xi,X2) ■b{xi)i'y^ c{x2), 

Jcc{x) — j j dXidX2 Fcc{x,Xi,X2) ■ c{xi)TccC{X2), 
Fbc{x, Xi,X2) = S{X - CiXi - C2X2)^bc{{xi - X2f) , 

Fcc{x, xi,X2) = S ^x - ^^^ ^ ^cc - 3:2)^) , 

Here wc adopt the notation: I — e~ , ^~ , t~ ,1 — e~^, , t'^ , ^ _ vrib / {mb+rric) C2 = mc/{mb+mc). 

The S-matrix element describing the semileptonic decays B^. {cc) + ID is written as 

Sb^^{cc) = 9bc9cc^Vbc JJJdxdXidX2B~{x)6{x-CiXi-C2X2)^bc{{xi-X2)'^) 



dydyidy2 6. 



Av) 5 [y - ^4^) *cc {{yi - 2/2)') j dz {10^ u)^ 



X ( T {b{xi) i 7^ c{x2) ■ c{y{) T,, c{y2) ■ c{z) O'^ c{z)] )o. 
The matrix element is calculated in the standard manner. We have 

TB^^(cc){Pl,P2,kl,K) = i{2TT)^S{p-i -p2-h- K) MB^^(cc){Pl,P2,kl,ku), 
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G 



1 

X -Tr 
4 



where pi and p2 are the -Be and (cc) momenta, respectively. The spin coupling structure of the (cc)-states is given by 

The calculation of the transition matrix elements proceeds along similar lines as in the case of the mass 
functions. For the scalar vertex functions one has 



She 



<^cc{-{k + P2/2f) 

Sg{k+p) = {mq+ ji+ ^) y"dae-"K-(fe+p)'). 



Again we shift the parameters ai{i = 1,2,3) according to 

ai {sbc + Scc)ai, 



I d^a f{a\,a2,as) = I dtt^ I d^a5 \l — ^^ai \ f{tai,tai,ta3) 



where d^a = daida2das. One then obtains 

M^{pi,P2) 
Ct 

Wbc 



1 r g(fe+cpj pi+cp3 P2)Vct i Xr(. • •) \ , 
J ■k'^ i 4 / ' 



1 



(sfec + Scc)(l+i)' 
C2 "^60 + tai 



1 + t 

Sbc 

■ 5 

Sbc y Sec 



^P2 



Wec/2 + ta2 
1 + t ' 



Occ 



where the symbol < ... > is related to the corresponding symbol < ... > defined in Sec. II Eq. (5). In the present case 
the symbol < ... > stands for the four- fold integral 



(•••) = (Sb, 



oo 1 / 3 \ 

+ scc)- Jdt-ij^ J d^asll-J2^n 



c)^ (. .. 



(•••), 



where 



z = t [{a2 + as) ml + ai - aictsPi - Q!2Q!3P2 - "i"2 q^] 

+ |pf [tWbcC2{2ai+a2-C2)+tWccOii/2-tai{ai+a2)+,WbcC2{wcc/'2.-C2+W},cC2)] 

+ pi, [tWbcC2a2 + tWcc{ai/2 + a2 - 1/4) -ta2 (ai + a2) + Wcc/4(2w6cC2 - 1 + Wcc)] 
+ [t{-'WbcC2a2-'WccCei/2 + aia2)-WbcWccC2/2] |. 
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One then shifts the loop momentum according to — ^- fc — Cp^ pi — Cp^ p2 ■ 

Our final results are given in terms of a set of invariant form factors defined by 

Be (cc) s=o ) = P'^ F+{q^) + {q^), (7) 

M>^{B,^{cc)s=i) = ^ et{-g^''PqAo{q') + P'^P-A+{q') + q^P-'A.{q^) 

+ie^''''f'Paq(3V{q^)}, (8) 

M^{Be^icc)s=2) = et„ {g''«p-Ti(g2)+p-p° [P>' T2{q^) + q^ mq^)] 

+ie^-''^P-Psq^T4{q^)}, (9) 

P = P1+P2, q = pi -P2- 

In our results we have dropped an overall phase factors which is irrelevant for the calculation of the decay widths. 

The calculation of traces and invariant integrations is done with help of FORM [27]. For the values of the model 
parameters (hadron sizes Ah and constituent quark masses niq) we use the values of [28]. The numerical evaluation 
of the form factors is done in FORTRAN. 

IV. ANGULAR DECAY DISTRIBUTIONS 

Consider the semileptonic decays B~{pi) — > {cc){p2) + l{ki) + i'{kv) and -B+(pi) — > {cc){p2) + l{ki) + i'{ki,). Recalling 
the expression for the matrix elements, one can write 

MB--.-cc{PuP2M,ku) = -^VbcM^{pi,p2)uf{ki)0^'v^;'iK), 



MB+^SciPUP2,kl,k^) = -^VbcMi,{pi,P2)ut{K)0^'vf{kl), 



where pi and p2 are the Be and (cc) momenta, respectively. 
The angular decay distribution reads 

_ G% ,2 (g'-/x')|P2| „ . . 

dg2dcos0 " (27r)3''^'=' ■ 8m?g2 ^^"^ 

where /U is the lepton mass and IP2I = A^/^(mf,m2,g^)/(2mi) is the momentum of the (cc)-meson in the Bc-rest 
frame. 

L^" is the lepton tensor given by 

L^^ = i Tr (O^ /^.O^ = fcf k-^ + k1 - 5"^ ±ie'^'''^f' ki^^k,^, (11) 

The lepton tensors L'^ and refer to the (iDi) and (li'i) cases. They differ in the sign of the parity-odd e-tensor 
contribution. The hadron tensor 7?^^ = Mfj,{pi,p2) A4l{pi,p2) is given by the corresponding tensor products of the 
transition matrix elements defined above. 

It is convenient to perform the Lorentz contractions in Eq. (10) with the help of helicity amplitudes as described in 
[29] and [30, 31]. First, we define an orthonormal and complete helicity basis e'^(m) with the three spin 1 components 
orthogonal to the momentum transfer g^, i.e. e^(m)g^ = for m = ±,0, and the spin (time)-component m = t 
with e^(t) = 

The orthonormality and completeness properties read 

el{m)e^{n) = Qmn (m,n = t,±,0), 

(12) 

<^n{m)el{n)gmn = 9ixv 
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with gmn = diag ( + , — , — , — ). Wc include the time component polarization vector e^{t) in the set because we 
want to include lepton mass effects in the following. 

Using the completeness property we rewrite the contraction of the lepton and hadron tensors in Eq. (10) according 

to 

= L{m,n)gmm'gnn'H{m',n') (13) 

where we have introduced the lepton and hadron tensors in the space of the helicity components 

L{m,n) = e''(m)e^"(n)V, if(m,n) = et''(m)e"(n)F^.. (14) 

The point is that the two tensors can be evaluated in two different Lorentz systems. The lepton tensors L{m,n) will 
be evaluated in the Iv or W-c.m. system whereas the hadron tensors H{m, n) will be evaluated in the Be rest system. 



A. Hadron tensor 



In the Br rest frame one has 



where 



V^t = (mi, 0, 0, 0), 

= 0, 0, -|p2|), (15) 

g'' = (90, 0, 0, |p2|), 



m\+m2- _m\-m\ + q^ 

-^2 = , <7o — 7, , 

2 mi 2 mi 

E2 + qo=mi, go=9^ + |P2|^ |P2p + ^^2 go = ^(m? - mi - g^). 



In the Bc-rest frame the polarization vectors of the effective current read 

e''(t) = ^(go, 0, 0, |p2|), 

vr 



e''(±) = -i=(0, Tl, -i, 0), (16) 
e'^(O) = ^(|P2|, 0, 0, go). 

vr 

Using this basis one can express the helicity components of the hadronic tensors through the invariant form factors 
defined in Eqs. (7-9). We treat the three spin cases in turn. 

(a) Be (cc) s=o transition: 

H{m,n) = (e+''(m)M,,) ■ (f^- („)>!„) ^ = HmHl (17) 
The helicity form factors Hm can be expressed in terms of the invariant form factors. One has 

Ht = -^{PqF+ + q^F-), 

vr 

H± = 0, (18) 

2mi |p2| 
^0 = 7^^+- 

Vr 
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(b) Be — > (cc) s=i transition: 

The nonvanishing helicity form factors are given by 

Hm = e^''{m)Mi,ael°'{m) for m = ±,0 (19) 

and 

Ht = e^^{t)M^o.el''{0) (20) 

As in Eq. (17) the hadronic tensor is given by H{m, n) = HmtlX- 

In order to express the hehcity form factors in terms of the invariant form factors Eq. (8) one needs to specify the 
heUcity components e2(m) (m = ±,0) of the polarization vector of the {cc)s=i state. They are given by 



e^(±) = -i=(0, ±1, -I, 0), 



£2(0) = — (|P2|, 0, 0, -£2). 
m2 



(21) 



They satisfy the orthonormality and completeness conditions: 



(22) 



The desired relations between the heUcity form factors and the invariant form factors are then 

= et^(t)e^"(0)Al^„ = 1 '^'\'^^ {p.q{-Ao + A+)+q^A_), 
mi +m2 mi^q^ 

H± = e+^(±)et"(±)M^„ = i {-P.qAoT2mi\p2\V), (23) 

mi + 7712 

Ho = e+^(0)et"(0)7W^„ 

^ ^ {-P-q{ml-ml-q^)Ao+4mi\p2\^A+). 



(c) -Be — * (cc) s=2 transition: 
The nonvanishing hehcity form factors can be calculated according to 

iJ„ = et^(m)Al^a,a,et"i"^(m) for m = ±,0 (24) 

and 

Ht = e^>^{t)M^o..ocJr'"\^) (25) 

Again the hadronic tensor is given by H[m,n) = H„^H^-^. 

For the further evaluation one needs to specify the helicity components e2(m) {m = ±2, ±1,0) of the polarization 
vector of the ccs=i- They are given by 

6r(±2) = 6^(±)6^(±), 

er (±1) = ^ (e^(±) £2(0) + e^(0) ,^^{±)) , (26) 
er(0) = -L(e^(+)e^(-) + e^(-)e^(+)) + y| 6^(0) 62^(0), 



11 



where ei^ir) are defined in Eq. (21). 

The relation between the helicity form factors Hm and the invariant form factors Eq. (9) read 

H, = h{t, 0) = ^|l + (|P2l' +E2qo + mi go) T2 + T3} , 

V O 1712 \/q 

H± = /i(±,±) = ^^^yPHl (riT2mi|p2|r4), 
v2m2 

H, = ft(0,0) = ./l^^^(K-m2-g2)Ti+4m?|p2pT2). 

V b V9 



B. Lepton tensor 

The hehcity components of the lepton tensors L{m, n) are evaluated in the {W)-c.m. system ki+kv = 0. One has 

= 0, 0, 0), 

k'i = (|ki|, |ki|sin6'cosx, |ki| sin6'sinx , |ki|cos6'), (27) 
ki = {El, — |ki| sin^cosx , — |ki| sin^sinx , — |ki|cos^), 



with 



In the (/P)-c.m. frame the longitudinal and time-component polarization vectors are given by 



e^{t) = ^ = (1,0,0,0), 



e''(±) = -L(0,T,-i,0), (28) 
e''(0) = (0,0,0,1). 

Using Eqs. (28) and (11) it is not difficult to evaluate the helicity representation L{m,n) of the lepton tensor. 

In this paper we are not interested in the azimuthal x distribution of the lepton pair. We therefore integrate over 
the azimuthal angle dependence of the lepton tensor. Of course, our formalism is general enough to allow for the 
inclusion of an azimuthal dependence if needed. After azimuthal integration the differential {q^,cos9) distribution 
reads 

-= - l + cos^e •— ^ + -sin2 0-— f TtCos^-— ^ 29 

dq^dcosO 8 ^ ' dq^ 4 dq^ A dq^ ^ ' 

where we take the polar angle to be the angle between the p2 and the ki in the lepton-neutrino cm. system. 
The upper and lower signs in front of the parity violating (p.v.) contribution refer to the two cases l~u and I'^u, 
respectively. 

The differential partial helicity rates dFi/dq^ and dTi/dq^ in Eq. (29) are defined by 

dfi _ ^2 
d^ ~ Vrfg^' 
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where we have introduced a standard set of hcUcity structure functions Hi {i = U, L, P, S, SL) given by the following 
linear combinations of the helicity components of the hadron tensor H(rn,n) = H^H^: 

Hu = Re (h+hI^ + Re [h_h1_^ Unpolarized — transverse 

Hl = Re (HoH^o) 

Longitudinal 

Hp = Re(^H+Hl^ - Re h1_ ) Parity - odd 

Hs = 3Re {HtH}^ Scalar 

HsL = Re (^HiHo^ Scalar — Longitudinal interference 

Note that the helicity amplitudes are real such that the complex conjugation symbol can in fact be dropped. 

It is evident that the "tilde" rates F in Eq. (30) do not contribute in the limit of vanishing lepton masses. In the 
present application this means that they can be neglected for the e- and /x-modes and only contribute to the r-modes. 

Integrating over cos one obtains the differential distribution 

dV _ dTu dTL dTu dTL dfs 
dq^ dq^ dq^ dq^ dq^ dq^ ' 

Finally, integrating over g-^, one obtains the total rate 

In Sec. V we list our predictions for the integrated partial helicity rates F^ (i ~ U, L, P) and Ti i ~ (U, L, S, SL). 

To save on notation in the following we shall sometimes use a self-explanatory notation for the differential and 
integrated partial helicity rates. For example, we write U for either the differential or the integrated helicity rates 
dTu/dq^ and Ff/, respectively, and U for dTu/dq^ and Tjj. 

An interesting quantity is the forward-backward asymmetry ApB oi the lepton in the [iD) cm system which is given 

by 

3 ±P + ASL 

Afb = - — 31 

+ U + L + L + S 

In Sec. V we shall give our numerical predictions for the asymmetry Aps for the decay channels under study. 

For the discovery channel Be J/xplv with the J/\p decaying into muon pairs the transverse/longitudinal composi- 
tion of the produced J/il) is of interest. The transverse/longitudinal composition can be determined by a measurement 
of the angular orientation of the back-to-back muon pairs in the J/V' rest frame. The relevant angular distribution 
reads 



dV 



dq^d cos 9* 




(32) 



where 6* is the polar angle of the muon pair relative to the original momentum direction of the J/ip. We have included 
lepton mass effects so that the angular decay distribution in Eq. (32) can be also used for the r-mode in this decay. 
The transverse and longitudinal contributions of the J/tp are given by {U + U) and {L + L + S), respectively. 

One can define an asymmetry parameter a* by rewriting Eq. (29) in terms of its cos^ 6* dependence, i.e. dV oc 
1 + a* cos^ 9* . The asymmetry parameter can be seen to be given by 

^.^f^+£-2(X + Z+|)_ (33) 
U + U + 2{L + L^- S) 

Our predictions for the asymmetry parameter a* appear in Sec. V. 

We have only written out single angle decay distributions in this paper. It is not difficult to write down joint 
angular decay distributions including also azimuthal correlations in our formalism if necessary. 
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V. NUMERICAL RESULTS 

Let us discuss the model parameters and their determination. Since we consider the decay of the Bc-meson into 
charnionium states only, the adjustable parameters arc the constituent masses of charm and bottom quarks and the 
size parameters of the B^-meson and charmonium states. The values of quark masses were determined in our previous 
studies (see, for example, [28]) of the leptonic and semileptonic decays of the low- lying pseudoscalar mesons (tt, K, 
D, Ds, B, Bg and Be). The values of the charm and bottom quarks were found to be to^ = 1.71 GeV and = 5.12 
GeV. The value of A^c = 1-96 GeV was determined from a fit to the world average of the leptonic decay constant 
fs^ = 360 MeV. The value of Aj/.^ = 2.62 GeV was found from a fit to the experimental value of the radiative decay 
constant = 405 McV which enters in the J /ij) e+e~ decay width (./^^J^* = 405 ± 17 MeV). 

In our calculation we are using free quark propagators with an effective constituent quark mass (see, Eq. 4). This 
imposes a very simple yet important constraint on the relations between the masses of the bound state and their 
constituents. One has to assume that the meson mass Mh is less than the sum of the masses of their constituents 

Mh < nig, + mq^ (34) 

in order to avoid the appearance of imaginary parts in physical amplitudes, which are described by the one-loop quark 
diagrams in our approach. This is satisfied for the low- lying pseudoscalar mesons tt, K, D, Dg, B, Bg, Be and 77c and 
also for the J/ip but is no longer true for the excited charmonium states considered here. We shall therefore employ 
identical masses for all excited charmonium states nice = mj/^=3.097 GeV (except for the ijc) in our matrix element 
calculations but use physical masses in the phase space calculation. This is quite a reliable approximation because the 
hyperfine splitting between the excited charmonium states and J/ijj is not large. For example, the maximum relative 
error is (m^(3836) - rnj/^)/mj/^ = 0.24. 

The size parameters of the excited charmonium states should be determined from a fit to the available experimental 
data for the two-photon and the radiative decays as was done for the J/ip-meson. However, the calculation of the 
matrix elements involving two photons will be very time consuming because one has to introduce the electromagnetic 
field into the nonlocal Lagrangians in Eqs. 1-3. This is done by using the path exponential (see, our recent papers 
[28] and [17]). The gauging of the nonlocal Lagrangian with spin 2 has not yet been done and is a project all of its 
own. For the time being, we are calculating the widths of the semileptonic decays Be — > (cc) -I- Iv by assuming an 
identical size parameters for all charmonium states A^c = Aj/^ = 2.62 GeV. 

In order to get a quantitative idea about the invariant form factors we list their q'^^^ = and g^^^ = ("^1 ~ ^^2)"^ 
values in Table II. 

We put our values of the decay rates in Table III together with those predicted in other papers. A number of 
calculations are devoted to the Be — » rjeliy and Be J/iplv decays. All of them predict values at the same order of 
magnitude. A study of the semileptonic decays of the Be-meson into excited charmonium states was done in [6, 21] 
within an approach which is quite different from our relativistic quark model. Concerning the electron-modes here is 
quite good agreement with [6, 21] in the case of Be J/^peu, rjcev, Xc2&i^- Our rates are a factor of 1.5 (1.8) larger 
for Be — * Xcoez^, heey decays and our rate is a factor 1.6 smaller for Be XciCJ^ decay. Concerning the r-modes here 
is quite good agreement with [6, 21] in the case of Be — > XcOTy, hcTV decays but our rates are almost a factor of two 
smaller for the other modes Be Xci^v, Xc2Tv. 

The partial rate for Be J/'ip -I- ^ -I- is the largest. The partial rates into the P-wave charmonium states are all 
of the same order of magnitude and are predicted to occur at ~ 10% of the most prominent decay Be ^ J/tp + 1 + u. 
The decays of the Be into D-wave charmonium state are suppressed. The r-modes are generally down by a factor of 
^ 10 compared to the e-modes except for the transitions Be rjc, Be J/if) and Be ■i/'(3836) where the r-modes 
are smaller only by a factor of ~ 3 — 4. 

In Table IV we list our results for the integrated partial helicity rates Fj (i = U, L, P, U, L, S, SL). They are needed 
for the calculation of the forward-backward asymmetry parameter Aps and, in the case of the decay Be J/ip + l + h', 
for the calculation of the asymmetry parameter a* determining the transverse/longitudinal composition of the J ftp in 
the decay. The partial "tilde" rates Fj are quite tiny for the e-mode as expected from Eq. (30) but are not negligible 
for the r-modes. This shows up in the calculated values for ApB in Table V. For the decays into spin states ApB 
is proportional to SL and thus tiny for the e-mode but nonnegligible for the r-modes. For the decay into the other 
spin states one has ApB{e~) = —ApB{e'^) but Apb{t^) ^ —Apb{t^) as can easily be appreciated by looking at 
Eq. (31). The forward-backward asymmetry can amount up to 40 %. The transverse and the longitudinal pieces of 
the J/ tp in the decay Be ^ J/ip + I + i' are almost equal for both the e- and the r-modes (see Table IV) . According 
to Eq. (33) this implies that the asymmetry parameter a* should be close to —33% as is indeed the case as the entries 
in Table V show. For the other two modes involving spin 1 charmonium states the transverse/longitudinal population 
is quite different. For the transition Be Xci the transverse mode dominates by a factor of ~ 3 for both the e- and 
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TABLE II: Predictions for the form factors of the Be — » (cc) tansitions. 





„2 
H 


-r 




'He 





0.61 


-0.32 




2 

'7m ax 


1.14 


-0.61 


XcQ 





0.40 


-1.00 




2 

?max 


0.65 


-1.63 









A+ 


A- 


V 







1.64 


0.54 


-0.95 


0.83 




2 

Qmax 


2.50 


0.97 


-1.76 


1.53 


Xcl 





-0.064 


-0.39 


1.52 


-1.18 




2 

Qmax 


0.46 


-0.50 


2.36 


-1.81 







0.44 


-1.08 


0.52 


0.25 




2 

Qmax 


0.54 


-1.80 


0.89 


0.365 







Ti 


r2,GcV-'' 


T3,GeV-2 


Ti, GeV"^ 






2 

'7m ax 


1.22 
1.69 


-0.011 
-0.018 


0.025 
0.040 


-0.021 
-0.033 






9max 


0.052 
0.35 


0.0071 
0.0090 


-0.036 
-0.052 


0.026 
0.038 



T-modes whereas for the transition Be he the longitudinal mode dominates by a factor of ^ 13 and ^ 7 for the e- 
and r-modes, respectively. 

Taking the central value of the CDF lifetime measurement t{Bc) ~ 0.46 • 10""'^^ s [1] and our predictions for the 
rates into the different charmonium states one finds branching fractions of ~ 2% and ^ 0.7% for the decays into the 
two S'-wave charmonium states J/tp and rjc, respectively, and branching fractions of ^ 0.2% for the deacys into the 
P-wave charmonium states. Considering the fact that there will be a yield of up to 10^° mesons per year at the 
Tevatron and LHC the semileptonic decays of the Be mesons into charmonium states studied in this paper offers a 
fascinating area of future research. 
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APPENDIX: CONVENTION FOR DIRAC 7-MATRICES AND THE ANTISYMMETRIC TENSOR IN 

MINKOWSKI SPACE 

We use the conventions of Bjorkcn-Drell. Thus we define the metric tensor and the totally antisymmetric e-tensor 
in Minkowski space by 5'"' = g^j, — diag(-|-, ) and £0123 = — e°^^^ = 1. For the partial and full contractions 

of a pair of £-tensors one finds 

'yfj.1yfj.2y ' fii'^ ti2y fJ3 ^ y fJiy fi^y IJ.3 



P ^viV2H3fJ.i _ -0(„''^ a"^ - a"^ a^^^ 

'^fllfl2fl3fJi'^ ~ ^ 5/(112 yfl-Lyfl2) 

p. c,Vxfl2H3fJi — a „Vi 

IJ.\li2H3lii'^ ~ ^yfxi 
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TABLE III: Semileptonic decay rates in units of 10"^^ GeV. We use \Vcb\ = 0.04. 





This model 


[6, 211 


[71 

L J 


[81 

L J 


[91 

L J 


[121 

L J 


[231 

L J 


R ^ T) PI/ 


10.7 


14.2 


11.1 


8.6 


11±1 


2.1 (6.9) 


5.9 


Be — > TjcT V 


3.52 






3.3±0.9 








R J/ih p ij 


28.2 


34.4 


30.2 


17.2 


28±5 


21 6 ("48 3"! 


17.7 


Be — .//' ( ■ r ;/ 


7.(S2 






7±2 










2.'>'2 


i.()(S() 












Be " \ rO 'i' '-^ 


l).2() 


0.21!) 












R — "v 1 p jy 


1.40 


2.206 












-Bf. — y Yri T z/ 


0.17 


0.346 












R _j. /j pi. 


4.42 


2.509 












Be hcTl/ 


0.38 


0.356 












Be -^Xc2ev 


2.92 


2.732 












Be ^ Xc^ru 


0.20 


0.422 












Be xp{383,6)eu 


0.13 














Be V(3836)ti' 


0.0031 















p. HllJ.2lJ.3H4 — _94 

<^ (1*1(1*2 /is (1*4'^ — 



We employ the following definition of the 7^-matrix 



5 = ^5 = i/7172^3 ^ ±e^,^,^3^,7'^i7''=7''3^^4 ^ / 



The Icptons with negative charge {I = e~, /i^, r~) are referred to as "leptons" whereas the positively charged 
leptons [ = e+, /Lt+, T"*" are referred to as "antileptons" . 
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TABLE IV: Partial helicity rates in units of 10"^^ GeV. 



Mode 


u 


U 


L 


L 


P 


S 


SL 


— > p v 








10.7 


0.34 10"® 





0.11 10"* 


0.34 10"® 


Be ^ rjcT v 








1.18 


0.27 





2.06 


0.42 


Be — ^ J/"^ 6 ^ 


14.0 


3.8 10"''' 
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0.31 10"® 


8.00 


0.88 10"® 


0.30 10"® 


Be —> J/lp T V 
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0.73 
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0.50 
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0.63 


0.31 


B^ — ^ v^n e u 

c /Leu ^ 
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0.11 10"® 





0.32 10"® 


0.11 10"® 


-Be — » Xco 
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TABLE V: Forward-backward asymmetry ApB and the asymmetry parameter a*. 



Mode 






a* 


Be ^ r)ceiy 


9.64 10"^ 


9.64 10"® 




Be If, T u 


().:',() 






Be -1 / 1 ' (' u 


().2i 


-0.21 


-0.34 


Be — .//(■-;/ 


0.29 


-I). 05 


-0.24 


Be — > XcO e v 


1.29 10"*^ 


1.29 10"'' 




Be XcO TV 


0.38 


0.38 




Be — > Xci e V 


-0.19 


0.19 




Be XclTV 


-0.24 


0.34 




Be heev 


0.036 


-0.036 




Be ^ he TV 


0.39 


0.30 




Be Xc2ev 


0.16 


-0.16 




Be Xc2TV 


0.32 


0.05 




Be V(3836)ei^ 


-0.21 


0.21 


-0.17 


Be ^(3836) T V 


-0.21 


0.41 


0.006 



